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Introduction 

Hawking radiation (HR), the spontaneous and steady emission of thermal radiation by black 
holes [lj, plays a crucial role in black hole thermodynamics [2H2] and in the attempts to construct 
theories of quantum gravity [6]. Given that these are still incomplete, it is of importance to 
determine under which conditions HR occurs, and how its properties depend or not on hypothesis 
concerning the ultra-violet behavior of the theory. In fact, because modes with arbitrary short 
wavelengths are involved in the standard derivation [7HS] , it was questioned [TU1 E] whether this 
process would still be present if Lorentz invariance were broken on very short scales. This question 
was motivated, on one hand, by the possibility that quantum gravitational effects could be effec- 
tively described by a non trivial dispersive relation, see e.g. II 2l. and, on the other hand, by the 
fact that phonons, or other collective degrees of freedom, obey dispersive equations at short dis- 
tance, e.g. when approaching the inter atomic scale. Hence dispersion must be taken into account 
when computing the phonon spectrum that a black hole analogue [7] would emit. To this end, 
Unruh [T3] wrote a dispersive wave equation in an acoustic black hole metric. He then numerically 
found that the thermal properties of the flux are robust, i.e. not significantly affected when A 3> k, 
where A is the ultra-violet dispersive frequency, and k the surface gravity of the black hole. In a 
subsequent numerical analysis [13], it was observed that "the radiation is astonishingly close to a 
perfect thermal spectrum" . This was confirmed to a higher accuracy in [15] , and partially explained 
by analytical treatments [T6H20] . 

In spite of these works, the origin of this astonishing robustness is not completely understood. 
This is due to our ignorance of the parameters governing the first deviations with respect to the 
standard thermal spectrum. In the present work, we show that when A 3> k the most relevant 
parameter is the extension of the near horizon region in which second order gradients can be 
neglected. 1 We reach these results by studying the validity limits of the connection formula [16- 
[19] encoding the scattering across the horizon. While the core of the calculation is based on 
the mode properties near the horizon - which are universal as they rely on a first order expansion 
around the horizon- the validity limits are essentially governed by the extension of the region where 
this expansion is valid. As a result, on the one hand, the leading order expressions are universal and 
agree with the standard relativistic ones, and on the other hand, the first deviations depend on this 
extension. Their evaluation is carried out using a superluminal dispersion relation. Interestingly, 
these results also apply to subluminal dispersion, as we show in establishing a correspondence 
between these two cases. We apply our analysis to both black and white holes, and show that HR 
in the latter backgrounds gives rise to undulations |24l [25] that contribute as classical waves to 
the observables. We then compare our predictions with the numerical works [26[ [27] which were 
done in parallel with the present analysis. We apply our treatment to the black hole laser effect in 

1 When these higher derivatives vanish, the background corresponds to a de Sitter space endowed with an ho- 
mogeneous preferred frame that respects some of the de Sitter isometries. This was noticed in [21) . and further 
developed with Jean Macher during his PhD |22| . This was exploited in |23| . and we hope to report on it soon. 
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App. [Bland to massive fields in App. [C] 

I. BLACK HOLE METRICS AND DISPERSIVE THEORIES 
A. The choices of settings 

To study the propagation of a dispersive field on a black hole geometry, basically two lines 
of thought can be followed. First one can study a particular fluid and derive the equations for 
linearized perturbations from the known microscopic theory, e.g. of a Bose condensate [28]. Sec- 
ondly, more abstractly, one can identify the ingredients that must be adopted in order to obtain a 
well-defined mode equation. We adopt here this second attitude [131 QUI [29] as it is more general, 
and as it reveals what are exactly the choices that should be made. 

First one needs to choose a dispersion relation in a flat (homogeneous) situation 

Q 2 = F 2 (p), (1) 

where fi is the frequency in the reference frame, p is the norm of the wave vector, and F 2 = 
p 2 db p 4 /A 2 + o(p 4 ). We thus work in units so that the group velocity is unity for O — > 0. The 
critical frequency A sets the scale at which dispersive effects become significant. In what follows, 
to simplify the equations we shall often work with 

^»=(p+|,) 2 =P 2 + £ + ^. 

At the end of the paper we shall consider more general relations. 

Second one needs to choose the two background fields, namely the metric g^ u , and the 'preferred' 
frame, i.e. the local frame in which eq. ([!]) is implemented. As explained in [2§] the latter can be 
described by a unit time-like vector field, it* 1 , in terms of which the preferred frequency is £1 = u^p^. 
In what follows, for reasons of simplicity, we work in 1+1 dimensions. Then the 'preferred' spatial 
wave vector p is p = s^p^, where s^u^ = 0, and s^s^ = —1. For further simplicity, we also work 
with stationary background fields, and we impose that the flow dx^/dr = is geodesic. Then, 
the Painleve-Gullstrand coordinates |30j . are the 'preferred' coordinates as they obey dt = u^dx^ 
and d x = s^d^. Using them, g^ u and are both encoded by a single function v(x): 

ds 2 = dt 2 - (dx -v(x)dt) 2 , (3) 

and 

17 = Lj — v(x)p, (4) 

where lj is the conserved frequency associated with the Killing field dt- Returning to the first 
attitude above discussed, both Eq. ^ and Eq. Q are obtained in the lab coordinates when 
considering the dispersive sound propagation characterized by Eq. ([l]) in a fluid whose flow velocity 
is v(x). 
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When v 2 (x) crosses 1, there is a Killing horizon. More precisely, for v < 0, it is a black 
(white) hole horizon when k = d x v evaluated at the horizon is positive (negative). This gradient 
determines the surface gravity of the horizon: neglecting dispersion, i.e. when F 2 = p 2 , the near 
horizon behavior of null right moving (upstream) geodesies is 

x = x e Kt . (5) 

Correspondingly, since xp = est near the horizon, the redshift experienced by right moving wave 
packets is [TBI 120] 

P = p e~ Kt . (6) 

The third choice concerns the quantization procedure. There exist indeed many inequivalent 
wave equations associated with the dispersion relation 

(u-vp) 2 = F 2 (p). (7) 

This remain true even when using an action to guarantee that the wave equation be governed by 
self-adjoint differential operators. In this paper we shall work with |13j 



[(u + id x v) (w + ivd x ) - F 2 (id x )] <p u = 0, (8) 

which is nothing but 

[(dt + d x v) (8 t + vd x ) + F 2 {id x )] <f> = 0, (9) 
applied on a stationary mode <fi = e~ v " t (p llJ . The associated conserved scalar product is 

{<h\<h) = i [ [fi(d t + vd x )(t> 2 - Mdt + vdjfa] dx. (10) 



Furthermore, we work with profiles v(x) defined on the entire real axis, and asymptotically constant: 
|u(±oo)| < oo. Hence the domain of integration is the entire real axis. The stationary (positive 
norm and real frequency) modes (f> u are then normalized by 

(<M<U = 6{(J-lS). (11) 

In conformity with second quantization [33], negative norm modes are named (0_ w )*, so that </>_ w 
are positive norm modes of negative frequency —oo. 

Eq. ^ reduces to the scalar massless d'Alembert equation when F 2 = p 2 = —d 2 . It differs in 
several respects from the (two dimensional) Bogoliubov-de Gennes equation [28]. It first differs in 
the hydrodynamical regime, i.e. in the limit A — > oo, in that the latter equation is not conformally 
invariant. As a result, left and right moving modes remain coupled in that case. Moreover, it 
also differs from Eq. ([8]) when taking into account the quartic dispersive effects. The differences 
arise from different orderings of d x s and v(x). To give another example, Eq. p]) also differs from 
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alternative models |16| [32] in which left and right moving modes remain exactly decoupled even 
when F{p) is non linear. 

Nevertheless, these wave equations share the same characteristics since these are determined by 
Eq. Q. What is less obvious is that these models also share, at leading order, the same deviations 
of the spectrum which are due to dispersion. This follows from the fact that these deviations 
are based on asymptotic expansions that are governed by Hamilton-Jacobi actions associated with 
Eq. @. 

B. in and out mode basis, connection formulae, and Hawking radiation 

In this paper, the properties of HR will be approximately determined by making use of con- 
nection formulae that relate asymptotic solutions of Eq. Q. Before describing this procedure in 
precise terms, let us briefly explain it. For the stationary profiles we consider, i.e. with v defined 
on the entire real axis and asymptotically constant, because of dispersion, the Bogoliubov trans- 
formation encoding the Hawking effect has the standard form of a scattering matrix. It should 
be stressed that this is not the case for relativistic fields. In that case indeed, because of Eq. <H, 
wave packets propagated backwards in time hug onto the horizon for arbitrary long time, and thus 
never transform into waves incoming from an asymptotic region. Instead, when there is dispersion, 
Eq. ([5]) is followed only for a finite time, and wave packets (propagated backwards in time) leave 
the near horizon region and reach, for superluminal dispersion, x = —00, see Fig. [T] 

Moreover when wave packets reach the asymptotic regions where v is constant, they can be 
decomposed in terms of stationary plane waves e ~ lu)t e % P^ x , Hence, the definition of the in and out 
modes is the standard one, see e.g. the scattering in a constant electric field |9]. The in modes 0™ 
are solutions of Eq. ^ such that the group velocity v gr = (ctap^) -1 of their asymptotic branches 
is oriented towards the horizon for one of them only. Hence when forming a wave packet of such 
modes, it initially describes a single packet traveling towards the horizon, whereas at later time 
it describes several packets moving away from the horizon. Similarly, the out modes (/>° ut contain 
only one asymptotic branch with a group velocity oriented away from the horizon. 

At this point it should be noticed that the dimensionality of these two sets depends on the 
asymptotic values of v. When v(x) contains one horizon, i.e. crosses —1 once, the dimensionality 
is 3 below a threshold frequency uj ma , x [E]: for w max > w > 0, there is one positive norm left 
moving mode (with respect to the fluid, but not the lab) 0^ ft , and a pair of right moving ones of 
opposite norm that we shall call and (^>_ CJ )* according to the sign of their norm. Hence the 
scattering matrix is 3 x 3. However, when v is smooth enough, 0j^ ft essentially decouples. This has 
been numerically shown in [15], and is mathematically justified in AppjAj Hence to a very good 
approximation, one recovers a 2 x 2 matrix characterizing right moving modes only. From now on 
we shall work within this approximation. 
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Figure 1: Space-time representation of the near horizon trajectories followed by wave packets made with 
right moving (with respect to the fluid) in modes in the presence of superluminal dispersion, e.g. that 
given in Eq. ([2]), and in the coordinate system (t,x) of Eq. We have also indicated their asymptotic 



amplitudes as given in Eq. (12 1. On the left panel we have represented the trajectories associated with 



and on the right those of the partner mode (4^ u )* with negative norm. The turning point of the latter is 



given in Eq. ( 28 ) 



Introducing the in and out sets of modes, they are related by 

(12) 





Because starred modes have a negative norm, the matrix is an element of U(l, 1). That is, the 
coefficients obey: 

I C^oj | | Pu) | — 1 7 

- Pu&uj = 0, (13) 

and |/3^| 2 = |Aj| 2 - Hence, when working in the in vacuum, the mean number of emitted pairs of 
quanta of opposite frequency is 

n u = \/3 u \ 2 . (14) 
In the relativistic limit, i.e. A — > oo, one gets the standard result 

-relativistic /-i c"\ 

U " ~ e 2nuj/ K _ ^ Li) > 

i.e., a Planck distribution at the Hawking temperature Th = n/2ir, in units where h = 1 and 
k B = I- 
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Our aim is to compute the coefficients of Eq. (12) using a connection formula. To this end, we 
first identify the various asymptotic solutions, and then we evaluate the globally defined solutions 
which we match to the asymptotic ones. These techniques have been already used in |16H19j . The 
novelty of our treatment concerns a careful control of the various approximations involved in this 



procedure. This will enable us to control the leading deviations from Eq. (15) which are due to 
dispersion, given v(x). 



C. The relevant properties of the profile v 

To be able to compute these leading deviations it is necessary to further discuss the properties 
of the profile v(x). When using relativistic fields, the temperature of HR is completely fixed by 
k = d x v, the gradient of v evaluated at the horizon, even though the asymptotic flux generally 
depends on other properties of v(x) which fix the 'grey body factors' [U[35]. However these describe 
an elastic scattering between </>^ ft and 4> u , and therefore do not affect the temperature, as can be 
verified by considering the equilibrium state described by the Hartle-Hawking vacuum (9j 134] . 
When dealing with dispersive fields, k is no longer the only relevant quantity. Indeed, as we shall 
show in the sequel, several properties of v{x) now become relevant. Moreover they govern different 
types of deviations with respect to the standard flux. For smooth profiles, there are basically two 
important properties, near horizon ones, and asymptotic ones. 

If there is a regular horizon at x = 0, v can be expanded as 

v{x) = -1 + kx + 0{x 2 ). (16) 

This near horizon behavior is only valid in a certain interval, not necessarily symmetric about 0. 
Hence we define D^ Q and D^ n such that for 

- <kx< LC (17) 

v is linear, to a good approximation (see region 1 in Fig. |2j). As we shall later establish, D^ n 
and control the leading deviations of the spectrum. It is worth noticing that in the limit 
D^ n ,D^ n — > oo, one effectively works in de Sitter space with a preferred frame since Eq. (|4j) and 
Eq. ([7]) still apply. In that limiting case, as we shall explain, the relativistic result of Eq. (15) is 
found with a higher accuracy. 

The other relevant parameter is related to the asymptotic values of v, that we assume to be 
finite. For superluminal dispersion, what matters is 

V (X = -OO) = -1 - < -1. (18) 

For subluminal instead, it is v(x = +oo) that matters. As discussed in [15] . Z) as determines the 



critical frequency w max (computed below in Eq. (29)) above which the flux vanishes 
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2 In that work, the profile was v = 



— 1 + D t&nh(Kx/D). Hence, when looking at the deviations of the spectrum 



s 



For finite values of x, v can have a complicated behavior. As we said, we only suppose that the 
interpolation between the asymptotic regions and around the horizon is smooth enough, so we can 
neglect non-adiabatic effects. Indeed, a sharpness in v(x) induces non-adiabatic effects |14 t 1151 [26] 
not related to the Hawking effect that both destroy the thermality of the spectrum and induce 
higher couplings between left and right moving modes, see Figs. 12 and 16 in |15j . These effects 
are due to the breakdown of the WKB approximation studied in AppfA} 



D. Hamilton-Jacobi actions and turning point 

1. Hamilton-Jacobi actions and turning point 

Eq. ([7]) can be interpreted as the Hamilton-Jacobi equation of a point particle. Indeed it suffices 
to consider the solution for p as a function of x at fixed u, which we call Pu,(x), as Pu,(x) = dxS^ 
where the action is decomposed as S = —ut + S u (x). One thus works the standard expression 

S u (z)= [ X Pu} (x')dx'- (19) 



In these classical terms, left and right moving solutions with respect to the fluid, i.e. the roots of 
u — vp = ±F, decouple and can be studied separately. Restricting attention to the right moving 



ones which enter in Eq. (12), one deals with 

u - v(x) Pu} (x) =F(p u (x)). (20) 
In the sequel, it will be useful to work in the p-representation with 

WUp) = -px + S u (x) = - J P X w {p')dp'. (21) 

In this representation, at fixed u), the position x is viewed as a function of p. It is given by 
X u = dpWu and obeys 

u-v{X UJ (p))p = F{p). (22) 

The usefulness of the ^-representation can be appreciated when considering the trajectories in near 
the horizon region, where v ~ — 1 + kx. In this region, irrespectively of F(p), one finds 

dp ( dX u 



at v a. ' (23) 

which gives Eq. Q, as in relativistic settings. Then the trajectory is algebraically given by 

x F (t) = X u (p(t)) = -- F{p) ~ P . (24) 
pn Kp 



upon changing D, the deviations associated with Dn n and D as were confused since both scaled in the same way. 
In fact, one of the main novelty of the present analysis, and the companion numerical works [261 127] , is to remove 
this confusion by analyzing the deviations due to D^ n and Dy ln only. 
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Unlike what is found for p(t), using Eq. ([2]), xf completely differs from Eq. ^ for \p\ > A. To 
adopt a language appropriate to the study of the modes, we shall work with to > only. Then 
negative frequency roots p- u > of Eq. (20) are replaced by the negative roots p u < associated 
with to > 0, as explained in [29]. Hence, Eq. (24) defines two trajectories, one with p > 0, and 
one with p < 0. At early times, i.e. for large \p\ 3> to and for superluminal dispersion \F/p\ > 1, 
both are coming from the supersonic region x < 0. Then, for p > the trajectory crosses the 
horizon and reaches x = +00, whereas for p < 0, it is reflected back to x = — 00, see FigJTJ What is 
important is that both trajectories stay in the near horizon region a finite time At. The integrated 
red-shift effect p final / Pinitial = e _KA * is thus finite, unlike what is found for relativistic propagation 
where Eq. ^ applies to arbitrary early times. 

For p < 0, the location of the turning point xt p is obtained by solving dx/dt = (d^p^ 
Using Eq. Eq. ^ gives 

n3 



0. 



to 



(1 + v) PuJ + 



PL 
2A 2 ' 



Hence, the momentum and the velocity at the turning point are (exactly) 

- Ptp = (10A 2 ) s , 



and 



v(x tp ) + 1 



(25) 



(26) 



(27) 



If to is sufficiently small, i.e. to < A (-Djf n ) 2 , the turning point is located in the near horizon region, 
and it is given by 

3 /ur 2 
2 \X 



KXtp 



(28) 



In classical terms, the main consequence of dispersion is the introduction of this turning point. It 



introduces a non trivial multiplicity of the real roots Pu(x), solutions of Eq. (25). This multiplicity 
will play a crucial role when solving the mode equation ([8]). From Eq. (27) and Eq. (18), we see 
that there is no turning point for to above 



a I Id, 



(29) 



This is threshold frequency w max mentioned in Sec IB It corresponds to the limiting case where 
the turning point x tp is sent to — 00. For to > to max only the positive root of Eq. (25) and the 
positive norm mode <$> w remain. Thus the transformation of Eq. (12) no longer exists. 

While these results have been obtained with a superluminal dispersion relation, however, they 



hold when the dispersion is subluminal. Indeed Eq. (25) is invariant under the three replacements: 

1 + -(1 + u), 

(30) 



to 
1 

A 2 



1 

"A 2 ' 
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1 + v(x) 




l.(b) l.(b) 



Figure 2: Shape of the typical velocity profile v, together with the extension of the relevant regions. For a 
given value of w, the near horizon region 1 where v ~ — 1 + kx splits into two: a region close to the turning 



point of Eq. (28) (l.(a)) where the WKB approximation fails, and two intermediate regions (l.(b)) where 
this approximation becomes reliable, and whose sizes are fixed by D(( n and D^ n . In the asymptotically flat 
regions 3, the solutions are superpositions of plane waves. The intermediate regions 2 play no significant 
role when v is sufficiently smooth, because the propagation is accurately WKB, i.e. no mode mixing. As 
we shall see, mode mixing essentially occurs at the scale of the turning point, i.e. in region l.(a). 

The first replacement exchanges the subsonic region and the supersonic one (for v < 0). As a result, 
a black hole horizon is replaced by a white one and vice versa. The second one amounts to a time 
reverse symmetry, t — > —t. At the classical level, it exchanges the roles of positive and negative 
roots of Eq. Q. At the mode level, it changes the sign of their norm, as discussed in Sec. IB 



The third replaces a superluminal quartic dispersion by a subluminal one. This exchange applies 
to all dispersion relations when expressed as F(p) — p — > —(F(p) — p). It replaces any dispersion 
that exhibits a superluminal character when p approaches A by the corresponding subluminal one. 
This correspondence thus applies to dispersion relations that pass from super to sub, as it is the 
case gravity waves in water when taking into account capillary waves [46J. 



Under Eq. (30), the trajectories are mapped into each other, as the function X w (p) of Eq. (22) 



is unchanged. Hence Eqs. ( 26|27p8 ), are also unaffected. Because we changed the sign of w, for 



subluminal dispersion, it is the trajectory of positive frequency that has a turning point. 
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2. The action So. 



th 



e near horizon region 



In preparation for the mode analysis we study the behavior of Su in the near horizon region 
where v is linear in x. Because of this linear character, it was appropriate to first solve the 
equations of motion in the p-representation and then go in the x-representation. This is also true 
for the action itself. Moreover, when solving Eq. Q, 4>ui{ x ) wm De obtained by inverse Fourier 
transforming the mode inp-space (f>u(p). Thus we express the action as = xp— W w (p). Imposing 
0, we get 



d p S 



Su(x,p ) = xp u (x) 



X UJ (p)dp, 



(31) 



I'D 



where p^ is a solution of Eq. ( 25 ) and po fixes the integration constant. Using Eq. ( 24 ) and Eq. Q , 
one gets 



Suj(x,p ) = xp L0 (x) ln(p w (x)) + 



6A 2 k 



+ e , 



where 9q is 



In (pq) 



Po 
6A 2 t 



(32) 



(33) 



To consider all solutions of Eq. ([8|) , we shall compute the action for all roots of Eq. ( |25| ) , including 
the complex ones. To this end, we need to define the integral dp' jp' = ]n(p/po), that arises 



from the first term of Eq. (24), for p complex. We shall work with the argument of cut equal to 
7r — e, so that ln(— 1) = — m. 



II. MODE ANALYSIS 
A. Asymptotic mode basis and x-WKB approximation 

Since the WKB approximation fails near a turning point, we cannot compute the coefficients 



of Eq. (12) using this approximation. In fact, under this approximation one would get a trivial 
result, namely f5 w = $ w = 0, |a w | = la^l = 1. To get a non-trivial result, we shall compute 
these coefficients by inverse Fourier transforming the modes in p-space and identifying the various 
terms sufficiently far away from the turning point, in a calculation that generalizes that of the Airy 
function (37]. 

In AppjA] we present the calculation of the WKB modes of Eq. ^ which generalizes the usual 
treatment since Eq. ^ is not second order. We also evaluate the errors made with respect to the 
exact solutions. In particular, the relevant errors are bounded by the inverse of the dimensionless 
parameter 

A(x) = A|2(l + „(x))|§. (34) 
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As expected, far away from the horizon A is large and the WKB approximation is accurate. More 



precisely A becomes of order 1 near x = xt p of Eq. (28) evaluated for to ~ k. Hence for these 
frequencies, which are the relevant ones for HR, A 3> 1 is reached for x/xt p S> 1. One also sees 
that at fixed x, A grows like A/re — )■ oo, hence A also governs the dispersion-less limit. 



The six roots p^ far away from the turning point 



Since we work in a weak dispersive regime, i.e. A> k, and since HR is related to frequencies 
to ~ k, we have to <C A for relevant frequencies. Moreover, since we impose that -D as is not too small, 
we also have to -C w raax , see Eq. (29). Hence u> max only concerns the high frequency properties of 



HR, which we no longer study. We focus instead on frequencies to ~ k. Even for such frequencies, 



the expressions of p u (x), solutions of Eq. (25), are complicated. However, their exact expression is 
not needed. It is sufficient to estimate them far away from the turning point, in order to build the 
mode basis. 

The denomination of the roots we use is based on that of the corresponding mode, which is 



itself based on the sign of the group velocity, as explained in Sec I B Moreover, we exploit the 
fact that for right moving modes the sign of the norm is that of the corresponding root p w (see 
e.g. [14 ). Therefore, for to > 0, positive roots shall be called p^, whereas negative roots shall be 
called in accord with the fact that negative norm modes are called see Fig. K3| Using 







■n 


PS * 




/ 












p 












Figure 3: Graphical resolution of Eq. ( 25 ) restricted to right moving modes. The dispersion relation F = 
p + p 3 1 2 A 2 and the straight lines fl = to — vp are plotted for a common value of to and two different values 



of v. On the left side, 1 + v < and we have 3 real roots given by Eqs. (35 36 37). On the right side, 



1 + v > and we have the single real root of Eq. (39). The two complex roots of Eqs. (40 41), are not 
represented. 
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this terminology, on the left of the turning point, one finds 



Pu 



P m , 



P°^ 



Ay/-2(1 + V) 

-Ay/-2(l + i 



00 



2(1 + v) 



00 



2(1 + 1;) 



:i + o(y)), 
(l + 0(y)), 



1 + v 



:i + 0(y 2 )), 



where the small parameter y is related to A of Eq. ( 34 ) by 

oj/k 
V= A(xj' 



(35) 
(36) 
(37) 

(38) 



Far away from the turning point of Eq. (28), y <C 1 and our expressions are reliable approximations. 
On the right (subsonic) region, one has only one real root 



ft 



out 



00 



1 + v 



(i + o(y 2 ))- 



(39) 



On this side there are also two complex solutions which do not correspond to any classical trajectory. 
However, when looking at the solutions of Eq. they govern the growing mode (f>1j and the 
decaying mode 4>w exactly as real roots govern WKB modes. These roots are 



pi = -iAV2(l + v) 
P i = iAy/2(l + v)- 



00 



2(1 + v 
2(1 + v) y 



,(i + o(y)), 

\ + 0{y)). 



(40) 
(41) 



On this side of the horizon, the corrections are again governed by y of Eq. (38). Hence, for to ~ k, 



the corrections to the roots are on both sides controlled by A of Eq. ( 34 ) . It should be also noticed 
that the errors for the two out roots are subdominant with respect to the other ones. 



2. The six WKB modes far away from the turning point 

In order to distinguish globally defined modes from their WKB approximations, the former 
shall be noted (f> u , and the latter ip^. Since we consider only 00 > 0, negative frequency modes (of 
negative norm) shall be written as (<p- w )* or (jp—J)*. 

Sufficiently far away from the turning point, i.e. for A S> 1, the WKB modes offer reliable 
solutions of Eq. (TsT) . They can thus be used as a basis to decompose globally defined solutions. 



We also assume that of Eq. (|17|) are large enough so that one can be at the same time be far 



away from the turning point and still in the region where v is linear. Using the expressions for the 



six roots in this region, Eq. (32), and neglecting the common phase depending on #0 of Eq. (33), 
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on the left side of the horizon one obtains 



vr- 



^ ~ [2A(x) 



[2A(a?) 



i e 



-i^ ln(AV2/s) 



^J\txk{\ + k\x\) 



e -^ln(Av^) 



/, nut\* 
IV-w) ~ e 



v / 4^tk(I+~k|x|) 

i i«— 



\/47ra; 



On the right one gets, 



out J^-i^ln(^) \ x 



7TW 



t 

<PL- ~ 



e ^ e -^ln(Av^) 

[2A(x)]3 - k\x\) 



e 2k e 



-i^ ln(AV2re) 



i i-i — --Adl 



ixr£ e §A(s). 



[2A(x)]5 V 4 ™^ 1 - K M) 
We now comment these expressions. Firstly, the relative errors of the two out modes are 



(42) 
(43) 
(44) 

(45) 
(46) 

(47) 



O 



A(x) 



whereas those of the four other modes are 



O 



1 + oj/k 
A(x) 



(48) 



(49) 



To get these expressions, we have taken into account two sources of errors: those coming from the 



approximate roots of Sec. II A 1 and those from the WKB approximation, see App. lA) They all 



depend on A of Eq. (34). Again, we see that the errors on the out modes, of low momentum p, are 
subdominant with respect to those of in modes which have a high momentum. As a result, Djj of 
Eq. (17) will be more relevant than . 



Secondly, whereas the normalization of the four oscillating modes is standard and based on the 



conserved scalar product of Eq. (10), those of the decaying ipz, and growing (pb modes follow from 



Eq. (A5) and the fact that of Eq. (32) is complex since the roots of Eqs. (40) and (41) are. 



Thirdly, since these two modes do not appear in the 'on-shell' Bogoliubov transformation of 



Eq. (12), we should explain why we are still considering them. First, the forthcoming connection 



formula will be a transfer matrix relating the general solution on each side of the horizon. Second, 
when considering problems with several horizons, these modes could contribute to the 'on-shell' 
S-matrix if they live in a finite (supersonic) size region between two horizons [36 . 
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B. Globally defined modes in the near horizon region 

The p-representation 

To accurately describe the behavior of the modes across the horizon one cannot used the x- 
WKB modes of the former Section. Rather, one should work in p-space, and look for solutions of 
the form: 

Mx)= [ 4>Up)e ipx ^=, (50) 

JC V27T 



where C is a contour in the complex p-plane. If it is well chosen, i. e. such that the integral converges 
and integrations by part can be performed, then it is sufficient that the dual mode <j> u satisfies 

(-u + pv) (-cj + vp) 4> u = F 2 (p)4> ul , (51) 



where v = v(x) = v(id p ). We should notice that Eq. (50) is a standard Fourier transform only if 
C is on the real line, something we shall not impose. The main interest of considering generalized 
contours is that it will enable us to compute all solutions of Eq. Q, including the growing ones. 

Because we only need the behavior of the modes near the horizon, the operator v in Eq. ( |51[ ) 
can be replaced by v = — 1 + indp. Hence one gets a second order ODE, irrespectively of F(p). 
The advantages to work in p-space are then clear \16\ [20] . Firstly, the solutions of Eq. ( |51| ) can be 
(exactly) written as a product 

M P ) = x(p)e-^ x^=, (52) 

V 47TK 

where \ obeys the w-independent equation: 

-K 2 p 2 d 2 pX = F 2 (p)x(p), (53) 

and where p~ % « _1 is independent of F. Hence, the deviations due to the dispersion F are entirely 
encoded in x- The origin of this factorization has to be found in the underlying structure of de 
Sitter space, see footnote [lj In addition, when considering the limit A — > oo in Eq. Q, p~ % ^~ x is 
exactly the relativistic (conformaly invariant) mode in p space. 



Secondly, unlike the original equation in x, Eq. (53) is perfectly regular. Moreover, when 
dispersion effects are weak, i.e. A//t > 1, the p-WKB approximation, so called not to confuse 
it with that used in the former Section, is very good. It generalizes what is done for the Airy 
function |37[ [30] where the mode equation in p-space is WKB exact. At this point, it is worth 



comparing the expression of the p-WKB modes in general and near the horizon. Using Eq. (21), 
one finds 



</w) ~ V—sT. n ^ ( 54 ) 



^^F(p) 



^4TTnpF(p) ^ VA 



l + Off)). (55) 
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The first line generalizes the standard WKB expression, and can be obtained by Fourier trans- 



forming Eq. (A5) at the saddle point approximation. It shows that (p w is universally governed by 
W u and X u , irrespectively of F(p) and v(x). The second line shows how exactly F(p) enters in 
<Pv in the near horizon region. In this region, since the mode equation is second order in d p , the 
corrections are uniformly boundedby k/A [40j. We also notice that using x* instead of \ in Eq. (52) 
would describe a left moving mode of negative norm [20 . This shows that the corrections to the 
p-WKB approximation describe creation of pairs of left and right moving modes, as in cosmol- 
ogy [15] . It is also of interest to notice that in models where left and right movers stays completely 
decoupled |16l 132] . the p-WKB modes of Eq. (55) are exact solutions in the near horizon region. 

We finally notice that Eq. (55), as the relativistic modep - *^ -1 , is well defined only when having 
chosen the branch cut of hip [T6HT9] . As explained below, different possibilities, and different 
contours C, lead to different modes. 



III. CONNECTION FORMULA 



A. The various modes in the near horizon region 

Using F{p) of Eq. ^ and Eq. (55) one gets 



1 f e (P - w ^> dp 

M&) = -nr= / 71 5Ti /s=- ( 56 ) 

V47TK Jc (! + ^2) 5 Pv2vr 

The forthcoming analysis generalizes former treatments for several reasons: 

1. Unlike |17H19j . we shall consider contours that are not homotope to the real line. This will 
allow us to obtain the general connection formula which includes the growing mode. 

2. We will make use of mathematical theorem of asymptotic development [40J under their exact 



form. This will lead to the proper identification of the validity conditions in Sec III C 



3. We will compute the phases of the Bogoliubov coefficients. These are essential to compute 
the correlation pattern (see Eq. (78)) and in the presence of several horizons, see AppjBj 



To evaluate Eq. (56), the first thing to take care of is the convergence of the integral. Indeed, 



. 3 

for large p, the dominant term is of the shape e ip , hence, we should impose that our contour C 
goes to infinity in regions where Im(p 3 ) > 0. The second step is to perform the integral at a saddle 
point approximation. For this, we make a change of variable such that the p and p 3 terms are of 
the same order for all values of A/k. We thus write p = Ay/2n\x\t and get: 

= 6 ^ \Z_ / 6 " .A W ( Slg n (a; )t + V) Jt (57) 
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Hence we see that the large parameter A(x) defined in Eq. (34) can be used to perform a saddle 



point approximation (see z in Eq. (58)). Therefore A will govern the deviations due to this 
approximation. 

For completeness we recall the saddle point theorem. If z is some large parameter: 

e iz f(Pj) 



A(p)e 



izf(p) 



dp 
2n 



y/-if"iPj)z 



1 + 



(58) 



where the pj are saddle points of /, i.e. f'(pj) = of smallest imaginary part, and where the 
square root takes its principal value. This formula is valid if and only if one can deform C such that 
Im(/(p) — f(pj)) is always positive [4Uj . The first correction Ej(f,A) involves higher derivatives 
of / and A evaluated at pj 

.5(/"') 2 A-3/"" f'A 



Ej(f,A) 



.A"_ J'" A' 

~ l J 7 + l W? 



12(/ 



//\3 



(59) 



1. The decaying mode 



We saw in Sec I D 1 that for negative frequency, the particle is reflected near the horizon, at the 



turning point of Eq. (24). Hence we expect that the corresponding mode will decay on the other 
side, for x > 0. This behavior is implemented by imposing that the branch cut is — iR+ and that 
the contour is homotope to the real line, as we now show. 



To evaluate Eq. (56) for x > 0, we use Eq. (57) and perform a saddle point calculation. The 
saddle points obey 1 + t 2 = 0. Just like in the Airy case [371 HO]) only t = i is relevant and its 
contribution is 



e 2«; 



g -i^ln(AV2K) 



2«e' 



•IA(«) 



l + O 



K\X\ 




(60) 



As required, the mode decays on this side of the horizon. The error term has been estimated using 



Eq. ( 58 ) . (More precisely, when computing Ej , we found a bounded function of oj/k and x times 



(1 + w 2 /k 2 )/A. This justifies our expression.) We notice that this expression coincides to (p^ of 



Eq. (46) up to a factor. Therefore, this factor defines the scattering coefficient and its x-dependent 



correction: 



l + O 




(61) 



We underline that this identification introduces no new errors because those due to the saddle 



point are of the same order as those already present in Eq. (46). Since the general method is now 



understood, we proceed with the same mode on the other side of the horizon, and then apply the 
same method for two other modes so as to get the general connexion formula. 

For x < 0, the saddle points now verify 1 — t 2 = 0. However, because of the branch cut, one 
must cut the contour into three separate branches, as shown in Figjlj C\ and C% go from ±oo + e 
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and dive toward — too on each side of the branch cut. C3 encircles the branch cut, and is necessary 
for the union of the 3 new contours to be homotope to the original one. Separating the three 
contributions, <p u = ip Cl + ip C2 + ip Cs , tp Cl and ip Ci are evaluated by the saddle point method and, 



after identification with the WKB modes of Eq. ( 42 ) and Eq. ( 43 ) , respectively give 



Ci 



vr- 



U1TT ■ 37T 

fe « e 4 



c 2 



1 + 




(62) 
(63) 



To properly evaluate (p C:i , one cannot use the saddle point method. However, because the factor 
e %px decays a l ng C3 , one can use a 'dominated convergence theorem', i.e. take the limit A — > oo 



in the integrand of Eq. (56). Using the Euler function, and </?°^ of Eq. (44) we get 

^2o7. 



As 




l + O 



— sinh ( — ) 

V k ) V 7r/ ; 

'«|a;|(l + ui 3 /k 3 



(64) 



The correction term has been calculated by expanding the integrand in Eq. (56) to first order in 
k/A and computing the integral again with the use of the T function. 




Figure 4: Representation of the contours in the complex p-plane determining the decaying mode, on the left 
pannel, when x is negative, and on the right when it is positive. The hatched regions are the asymptotically 
forbidden ones, and the dots indicate the saddle points that contribute to the integral. 



2. The transmitted mode 



To get another mode, we keep the same contour but the branch cut is now taken to be As 
we shall see it corresponds to a transmitted mode. For x < the saddle points still obey 1 — t 2 = 0, 
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but we can now use the saddle point approximation because the branch cut is no longer in the way. 
Taking into account that on the negative real t-axis lni = In \t\ — itr, we get: 



x e*4 + 



<P- 



UJ7V ■ 37T 

x e » e 4 



1 + 




(65) 



On the other side, for x > 0, the saddle points obey 1 + t 1 = 0, as previously. Because of the 
branch cut one cannot pick the contribution of the decaying saddle point. We must instead deform 
the contour to a region where the 'dominated convergence theorem' can be used and then stick it 



to the branch cut. With a computation similar to what was done for Eq. (64), we find 




U) 



e 2« e 



(66) 



3. The growing mode 



To get a third linearly independent solution of Eq. Q, we must construct the growing mode. 
To get it, we re-use the above defined contours C\ and C2, and we choose the branch cut to be 



For x < 0, the relevant saddle points are —1 for C\ and +1 for C2, and respectively give Eq. (62) 



and Eq. (C8). For x > instead, for both contours, the relevant saddle point is t = —i, and this 



gives 



P. 



£2 



-e « (ft 1 



l + O 



1 + 



73 



A(x) 



(67) 



Since all combinations of ip 1 ^ 1 and <^£ 2 behave as iplj when A — > 00, there is an ambiguity in choosing 
the third mode we shall use. We appeal to the conserved Wronskian to fix the choice. For a forth 
order differential equation, the Wronskian is a 4 x 4 determinant, but because we neglected the 
v-modes, it becomes 3x3. Once we have chosen 2 propagating modes, there is a unique choice of 
growing mode such that the basis has a unit Wronskian. The connection matrix is then an element 
of the group SL^(C), i.e. of unit determinant (up to an overall gauge phase). 

Therefore, our third mode is ip^j = i^£ 2 — e «" ' p^J . On the right and on the left we respectively 
get: 



x>0 



x<0 



1 + 




LUTV 

' « e 



1+0 



i+Jz 

A(x) 



(68) 
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B. Connection matrix and on-shell Bogoliubov transformation 



1. The connection formula 

The results of the former subsection can be synthesized by the following 3x3 'off-shell transfer 
matrix' that connects the WKB modes defined on either side of the horizon 



\<pt/ 



( 



(U Bn ) 



\ 



(69) 



{A 

We define C/bh through its transpose so that it relates the three amplitudes of any mode decomposed 



on the left and right side basis of Sec II A 2 Ignoring for the moment the correction terms, the 
matrix is 



Ubu 



e*5 \ 
2 



e * e 2] 



: 7T UJ-K 



\ 







— e 4 e « 



e k 



e 



(70) 



To simplify the above, we defined the 'normalized' T function: 

t(z) = T(-iz 

which obeys for large z 



sinh(7rz)e-^e i2ln(2) - i2 e-^, 



1 - e 



-2-kz 



12z 



T(z) 
Arg (t(z) 

As expected from our choice of modes, the determinant is a pure phase: 

det(C7 BH ) = e^- 



(71) 

(72) 
(73) 

(74) 



2. Robustness of black hole radiation 



Using Eq. (70) we can now easily extract the Bogoliubov coefficients of Eq. (12). Let us start 
with (f>^. Being a physical mode, it is asymptotically bounded, and therefore the amplitude mul- 



tiplying the growing mode should vanish. Moreover, on the left it asymptotes to (p™ of Eq. (42). 
Hence its six amplitudes obey 







Ubh 



x<0 



(a \ 



(75) 



z>0 
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where d u is the amplitude of the decaying mode. From these equations, and the corresponding 
ones for the negative frequency mode 



, the coefficients of Eq. (12) are 



— e « 
5. 



The amplitudes of the decaying modes in </>™ and 



• 2 



(76) 



are also fixed and given by 



2sinh(^f)' 
1 

2sinh(^f)' 



(77) 



For uj <C w max of Eq. (29), and when ignoring the x-dependent corrections of the former subsec- 



tion (i.e. to leading order in ft/ A), the mean occupation number Eq. (14) is exactly the relativistic 



result of Eq. (15). This is in agreement with what was found in |16H19| . although the conditions 
are now stated more precisely. This result implies that the spectral deviations due to dispersion 
are to be found by examining the various approximations that have been used. 

In addition, it should be noticed that sufficiently far away from the horizon, i.e. in a region 



where A(x) of Eq. (34) is larger than 1, the space-time correlation pattern of the Hawking particles 



of positive frequency and their inside partners of negative frequency is also unaffected by dispersion. 
This second aspect of the robustness of HR can be established by either forming wave packets of 
in-modes [16], i.e. considering non- vacuum states described by coherent states, see App.C. in [28J, 
or by computing the 2-point correlation function G(t,x; t'x') = {4>(t, x) <p(t' , x')) [4"lTl4"3"] . For a 
comparison of the two approaches, see |23j . In the in- vacuum, at equal time, the oj contribution of 
G for x > and x' < is given by 



G u [x, x') = a^l y™\x)y"t(x>) + u a* u {^(x))* (¥>-^')) S 



(78) 



see the B w term in Sec. IV. F. in |28j . Using the expressions of Sec II A 2| together with Eq. (76), far 

but still 

~Rie\x/x' 



away from the turning point of Eq. (28) but still in the near horizon region, we get 

'1 +oj 2 /k 2 



Guj(x,x') 



1 



1 + 



(79) 



sinh(^) 4vrw ' V~ ' V A 0*0 

At leading order in k/A, this is exactly the relativistic result. Hence, the long distance correlations 
are also robust when introducing short distance dispersion. This follows from the fact that the 
phase of /3 w /a w (and not just its norm) is not modified by dispersion. At this point we should 
say that this phase actually depends on those of the in and out modes that can be arbitrarily 
chosen. Thus, as such it is not an observable quantity. However, the phase of a^/3^(/j° u V°^ is an 
observable, independent of these choices. We have chosen to work with in and out bases where 
all modes have a common phase at a given p, see Sec ID 2 as this ensures that arg(/3 w /a w ) is 
unaffected by this arbitrary phase. 
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Using these basis, the phases of a u and of Eq. (76) also have a clear meaning as can be 
seen by considering the scattering of classical waves. They characterize the phase shifts which are 



not taken into account by the WKB modes of Sec. II A 2 In fact, using Eq. (71) one verifies that 
in the limit uj/k —> oo one recovers the standard WKB results, i.e. arg(a w ) = arg(r(o;/«;)) —> 
for the transmitted mode, and arg(a aJ ) — )■ tt/2 for the reflected one. For smaller values of oj/k, 



&rg(T(uj/ k,)) thus accounts for the non-trivial phase shift. In App. [Bjwe show that in the presence 
of two horizons, this shift affects the spectrum of trapped modes. 



3. White holes 

It is rather easy to apply our results to white holes as they behave as the time reverse of black 
holes. In terms of stationary modes, the correspondence is made by exchanging the role of in and 
out and perform a complex conjugation: 



-in.WH 



-out.BH 



(80) 



as in Eq. (D3) of [28]. This follows from a symmetry of Eq. Q. Indeed, the mode equation is 
invariant under 



-u) and v 



-v. 



For more details, we refer to the App.D of [28]. Here, it implies 

Uwr = (^bh)* , 



(81) 



(82) 



where [%h is defined through the same equation as Eq. (69). This implies that Bogoliubov 



coefficients of a white hole posses the same norm as those in the corresponding black hole setup. 

It is more interesting to look at correlations. Indeed, since both out modes now live on the same 
side of the horizon and contain high momenta, the correlation pattern drastically differs from the 
black hole one |24j . To describe this pattern, in order not to introduce confusing notations, the 
modes shall be called according to their status in a black hole geometry, that is, in the following 



expressions the modes are those of Sec II A 2 Then, in the white hole in- vacuum, the w-contribution 
of the equal time correlation function reads 



G™ H (x, x') = 0° ut (*) (CV))* + <TM («*'))* • 
Decomposing the BH out modes in their in content, the first term gives 



(83) 



(CV))* 



|2 CW (CV))* - otJMx)4^{x% 
K (CJ*))* (CV))* + l/il 2 {fatoY fc>V). (84) 



a. 



The most interesting phenomenon occurs when one looks at the zero frequency limit. In this 
regime, on the one hand, the zero frequency limit of (f>™ and (p 1 ^^ agree with each other, as can be 



23 



seen in Eq. (42) and Eq. (43), and equal o n E C which is a non-trivial function of x. On the other 
hand, when to 



0, Eq. (|76|) gives 

2 



2ituj 



and 



2ituj 



Then, since both terms of Eq. (83) agree in the limit, we finally get 

4k. 



G 



WH 



(x,x') ~ — 

m0 



Re ie^l <j$(x)\ Re {e"^ $V)} 



(85) 



(86) 



Before commenting on this expression, it is interesting to compare it with the corresponding one 
obtained in a black hole geometry. In that case, when working in the black hole in vacuum, the to 



contribution of the two point function is, see Eq. (83), 

gT = C(*) (<€(*'))* + (<P^')Y 

Hence in the limit to — > it gives (see Eq. 32 in 



G, 



BH 



2<(x) (4V))* 



(87) 



(88) 



which behaves very differently than G WH of Eq. (86). It does not diverge as 1/to and it is not the 
product of two real waves. To get an expression that might correspond to that of Eq. (86) one 



should express the in modes in terms of the out black hole modes given in Eq. ( 44 ) and ( 45 ) , as 
done in Eq. (78). Doing so, the prefactor of Eq. (86) is recovered but the spatial behavior of G BH 



is completely different because the out modes of Eq. (44) and (45) are defined on opposite sides of 
the horizon, and become constant in their domain. Hence unlike G WH , G BH cannot be written in 
the limit to — > as a product of twice the same real wave. 

In Eq. (86), we see that C? WH factorizes in the zero frequency limit as the 2-point function 
in inflationary cosmology when neglecting the decaying mode, see e.g. Eq. (20) in 07]. This 
fact shows that Re |e _ *4 0g n (a;)| contributes to G WH as in a stochastic ensemble of classical 
waves, i.e. each member of the ensemble contributes like a coherent state with a given phase, 
see App.C in [28J, but the amplitude of 

primordial fluctuations in cosmology. Moreover, the real function Re |e 



q 1 is still a Gaussian random variable, again as for 



^x) j corresponds to 

the undulation observed in [25], and the present analysis shows that this phenomenon is directly 
related to the Hawking effect. However, what fixes the constant amplitude found in the experiment 
remains to be understood because, when summing Eq. (86) over to, the factor 1/oj engenders a 
logarithmic growth in ln(t) as in Ref. |24j . Eq. (86) differs from the expression of [23] (which also 
factorizes) in that the surface gravity has not been sent to oo. Using Sec II A 2 one can compute 



its profile (in the WKB approximation) in the regions l.(b) of Figj2j In the right region l.(b) it 



decays according to the zero frequency limit of Eq. ( 60 ) , whereas in the left region, one has 

I 



Re 



(*)} 



cos(§A(x) + f ) 



(89) 



V&KK y / A(x)(l + k\x\) ' 

On the right it thus behaves very much like the decaying Airy function Ai \67\ 14*0] , and on the left 
it oscillates in a similar manner but, quite surprisingly, the phase shift 7r/4 has the opposite sign. 
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The origin of this flip is to be found in the extra factor of 1/p in the integrand of Eq. (56) that 
is associated with the relativistic (non-positive) norm of Eq. (10). It would be very interesting to 
observe the profile of Eq. (89) and its unusual phase shift in future experiments. Further away 



from the horizon, the undulation profile can be obtained from the zero frequency limit of Eq. (A5). 



C. Validity of the connection formula 

Our computation is based on two approximations. The first one is the p-WKB approximation 



introduced when solving Eq. (53) in the near horizon region. Its validity requires 

- » 1. (90) 

K 



This condition is the expected one. It involves neither uj nor the parameters Dy m of Eq. (17). 
Moreover, as we shall see below, it will not be the most relevant one in the general case. This 
is a non-trivial result. In addition, since the corrections to this approximation mix left and right 
moving modes [20], at leading order, the same spectral deviations will be obtained when considering 
models [HI [32] where the decoupling between these modes is exact. It would be interesting to 
validate this prediction by numerical analysis. 



The second approximations are controlled by A of Eq. (34). This quantity governs both the 



validity of the saddle point approximation, as in Eq. (60), and that of the WKB modes of Eq. 



(42 47). Since these corrections decrease when A increases, the pasting of the near horizon modes 
on the WKB ones should be done at the edges of the near horizon region. One could imagine 
pasting the modes further away, but this would require to control (j>w(p) outside the region where 
v is linear in x, i.e. to deal with ODE in d p of order higher than 2. This is perhaps possible but it 
requires other techniques than those we used, see [50] for recent developments. In any case, what 
we do is sufficient to control the error terms in the relativistic limit, and more precisely to find an 
upper bound. 

Being confined to stay within the near horizon region, the validity of the pasting procedure 
requires that 

A p = A(x past ing) = -(KX pastin g)5 ~ — (£> lin )2 > l ) (91) 

where D\\ n characterizes the extension of this region on either side of the horizon. As we see in 
FigjSJ at fixed A/re, the spectrum is very close to the relativistic one of Eq. (15) if D\[ n is large 



enough. Instead, below a certain threshold, the deviations become non negligible. Our criterion in 



Eq. (91) indicates that this will happen when A p is of order 1. Hence the threshold value for Dn r 



should scale as (re/A) 2 / 3 . This prediction is confirmed by the numerical analysis of |26| . 



We can be more precise. Indeed, as discussed after Eq. (16), the near horizon region is not 
necessarily symmetric. Hence, the values of ^pasting on the right and on the left will be different. 
This is important because error terms coming from in modes are dominant compared to those from 
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Figure 5: Deviation of the temperature found using the code of |26j for various flows v(x) which all have the 
same surface gravity k. Th = n/2ir is the usual Hawking temperature. To is the actual temperature given 
quartic superluminal dispersion computed for uj <C k. The parameter k p (see [26] for its precise definition) 
characterizes the slope of v(x) outside the near horizon region, i.e. in regions 2 of Figj2j At fixed A/n = 15, 
the critical value of Dn n ~ 0.2 below which Tq deviates from Th does not depend on k„. This is in agreement 



with Eq. (91) 



(92) 



out modes, see Sec II A Since the former lives on the left side, the validity condition is 

aj~£(i&) s >i, 

hj 

The higher sensitivity of the spectrum to perturbations of v localized on the left hand side was 
clearly observed in [26], see Fig. 8 right panel. This sensitivity has been recently exploited in [Hj 
to produce resonant effects. We can now estimate the deviations on the spectrum. Considering 
Eqs. Q48|49|60|64|92"l ), we obtain 
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where P is a polynomial function of degree 2. Therefore, at fixed uj < k, i.e. in the relevant 
regime for Hawking radiation, the leading deviations are bounded by a quantity scaling as 1/A^. 



In particular, Eq. (93) demonstrates that the spectrum is thermal for arbitrary small frequencies, 
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contrary to what has been claimed in [17\ 131] . Moreover, our analysis indicates that the deviations 
should grow with oj. This is compatible with the fact that n vanishes for oj > w TOax defined in 
Eq. (29), irrespectively of the value of the ratio A/k. However, we have not yet been able to 
confirm this precisely with the code of |26j . One of the reason is that the spectrum is radically 
modified when u approaches w max - So far the p-WKB approximation has been a subdominant 



effect. However, if one takes the de Sitter limit, i.e. Dn n — > oo, the correction term in Eq. (91) 
vanishes and p-WKB becomes the only source of deviations. 

To conclude we recall that the deviations have been computed using Eq. ([8]). In other mode 
equations, like that for phonons in a BEC |28j . the corrections to the x-WKB approximation will be 
in general larger. However these corrections are not due to dispersion but rather to the fact that the 
conformal invariance of Eq. ^ in the dispersionless limit will be lost. These corrections can thus be 
studied without introducing dispersion. This is also true when introducing an infrared modification 
of Eq. ^ associated with a mass or with a non- vanishing perpendicular momentum. In AppjCj 
we show that the spectral properties are still robust in that Eq. (92) is sufficient to guarantee that 
to leading order the Bogoliubov coefficients are unaffected by short distance dispersion. 



D. General superluminal dispersion 



Instead of Eq. ((2J), we now consider 

F\p) 



P + 
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,2n+l \ 2 



A 2r 



(94) 



taken again, for simplicity, to be a perfect square. In Sec II B we computed p-WKB modes for any 



dispersion in Eq. (|55|), hence the globally defined modes for F satisfying Eq. (94) are, see Sec. Ill 

dp 
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There now exist In + 1 linearly independent modes. Hence in terms of contours, there are 2n + 1 
sectors (Stokes lines) toward oo in the complex p plane. By using a contour homotope to the real 
line, and the same two possible branch cuts of lnp on ±iIR + , we can compute the 2 'on-shell' modes 
that are asymptotically bounded. Even though, there exist n pairs of growing and decaying modes 
on the subsonic side, and n — 1 pairs on the other side, only one pair in the subsonic sector will be 



relevant in the 'off-shell' connection formula of Eq. (70). Indeed, all the others pairs do not mix 



with propagating modes. Therefore, the different contours giving rise to the relevant modes will 



be quite similar to those of Sec III 



To perform a saddle point approximation, we introduce t = p/A\nx\ 2™ and get: 
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By a computation similar to that of Sec. Ill, at leading order in k/A, we recover the Bogoliubov 



coefficients of Eq. (76), thereby establishing their robustness for arbitrary integer values of n. 



Moreover, the deviations from this result are now governed by 

A T 2n + 1 

rv 

Hence, the error on the mean number of emitted quanta satisfies 

;relativistic / 
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E. Relating subluminal dispersion relations to superluminal ones 



So far we analyzed only superluminal dispersion relations. We should thus inquire how would 



subluminal dispersions affect the spectrum. At the classical level, as noticed at the end of Sec. I D 1 



there is a exact correspondence between these two cases. At the level of the modes Eq. (30) does not 
leave Eq. Q invariant as it does not apply to the left moving solutions governed by oj — vp = —F. 
However, it becomes a symmetry when neglecting the mode mixing between left and right movers. 



Therefore, in models where the decoupling between these is exact [161 E2]> Eq. (30) is an exact 
symmetry. Moreover, since the mode mixing between left and right movers is subdominant for 
general mode equations, the discrepancy of the spectral deviations between superluminal and 
subluminal dispersion will not show up at leading order. This is precisely what has been observed 
in Sec.VI.2 of [15]. 



At the level of the connection formula of Eq. (70), the three exchanges of Eq. (30) still are 



an exact symmetry since the [TeH-matrix is based on the right moving mode of Eq. (55) which is 



determined by the action W u (p). Therefore, at leading order in k/A, without any further calcula- 
tion, this symmetry implies that the spectrum of HR is equally robust for subluminal dispersion. 
Moreover, the leading order deviations from the thermal spectrum will be governed by the same 



expression as Eq. (98). It should be noticed that when applying Eq. (30), D^ n characterizes the 
extension on the sub-sonic region. 



The above symmetry should not be confused with that of Eq. (81 ) which also relates black and 



white holes. Indeed, the latter exchanges the role of left and right moving modes, while the former 
applies only to the right moving sector. Instead, these two symmetries can be composed with each 
other. This allows to compare black hole spectra without referring to white holes. 

To conclude the discussion, we mention that this approximate symmetry allows not only to 
predict several effects, but also to predict how the observables will quantitatively behave: 

• a laser effect will be found for sub-luminal dispersion in a flow possessing two horizons that 
passes from super to sub and then back to a supersonic flow, 3 and this exactly for the same 



3 We are grateful to Daniele Faccio and William Unruh for bringing our attention to this possibility in a discussion 
that took place during the Nice Colloquium on 'Analog gravity' in June 2010. 
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reasons that the laser effect was found in the 'reversed' flow in the case of superluminal 
dispersion [51] and AppjB} 

the frequencies and the growth rates of this subsonic laser effect will be governed by the same 
expressions as those of AppjB] (when neglecting the coupling to the left moving modes). 

subsonic phonon propagation in a non- homogeneous flow that remains everywhere subsonic, 
i.e. without a sonic horizon, will be governed by a 4 x 4 S-matrix that encodes new pair 
creation channels with respect to those found in the presence of a sonic horizon [50] exactly 
for the same reasons that a supersonic phonon propagation in a non-homogeneous flow that 
remains everywhere supersonic does so, as mentioned in [27] . 

the behavior of the Bogoliubov coefficients in the two cases will behave quantitatively in the 
same way. 

the undulation observed in white hole flow for subluminal gravity waves in the experiment 
of [25] is generated for the same reasons as that found in white holes for Bose condensates 
using the Bogoliubov-de Gennes equation [53] . As we discussed above, this is obtained by 



composing the symmetries of Eq. (81) and Eq. (30). 



Moreover, when the dispersion relations and the profiles v{x) + c(x) (where c(x) is the speed 



of sound that generalizes the 1 in 1 + v) obey Eq. (30) up to a possible rescaling of A, the 



momentum p, and distances, these undulations should have the same spatial profile. 



IV. CONCLUSION 

In this paper, we described the scattering of a dispersive field on a stationary black hole horizon. 
We computed the connection formula which relates WKB modes on each side of the horizon 
sufficiently far away from the turning point. Our main results are: 



Eq. (70) applies 'off-shell', which means that the contribution of the growing mode is taken 



into account. When requiring that its amplitude vanishes, Eq. (70) fixes both the 'on-shell' 



2x2 matrix of Eq. (12) encoding the Bogoliubov transformation between physical modes, 



and the amplitude of the decaying mode. In situations with several horizons, the presence 



of a growing mode in a finite size region could alter the scattering. In that case, Eq. (70) 
should be used to estimate this effect. 

The phases of the scattering coefficients are computed. This allows us to show that not 
only the out-going flux is robust against dispersion, but also the correlation pattern between 



Hawking quanta and their partner (in Eq. (76) and Eq. (78)). Moreover, in situations with 



several horizons, such as that of the black hole laser, these shifts enter in interference effects 
and directly affect the observables. Our prediction for their values is numerically confirmed 
as discussed in AppjB] 
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We characterized the leading spectral deviations due to dispersion. These critically depend 
on the spatial extension of the near horizon geometry, and not only on the ratio k/A. Indeed, 



they are governed by the parameter D^ n of Eq. ( 17 ) , and this for for a large class of dispersion, 



see Sec III D These spectral deviations imply only an upper bound on the frequency, in 
contradiction to what was claimed in \17\ [3Tj . 

A symmetry between superluminal and subluminal dispersion relations establishes that the 
spectral deviations due to one type of dispersion will be also found for the corresponding 



one when modifying the flow v(x) in a certain manner, see Eq. (30). 



In App. [Cj we extend our results to massive fields. In that case, the dispersion relation 
is modified both in the infrared and the ultraviolet. This induces two types of deviations 



from Eq. ( 76 ) . We show that these deviations decouple when m <C A and can be bounded 



separately. The ultraviolet effects are still controlled by Eq. (91) while the infrared ones can 
be studied with the massive relativistic equation. 
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Appendix A: WKB dispersive modes 



1. The method 



To apply a WKB approximation to the solutions of Eq. ^ , we write the mode as 
Injecting this into Eq. (|8|), we obtain 

"1 



(Al) 



(UJ — v\x) 



) 2 -F 2 



-id x 



-d k F z {K 



Ak^d^k^ + <i{d x k UJ ) 2 id x k u 
A2 



(A2) 



For definiteness and simplicity, the last term is given for F(k) = k 2 + fc 4 /A 2 but can be generalized 
to any polynomial dispersion relation. So far, this equation is exact. It is known as a Riccati 
equation [40] and was already used in the present context in [17J. It is adapted to a perturbative 
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resolution where the different terms are sorted in order of derivatives 
Hence we write kn, as 



k _ 1.(0) , k (i) , U2) 



+ 



here spatial gradients. 



(A3) 



where the superscript gives the number of derivatives (one way to sort the terms of Eq. ( A2 ) is to 
make the scale change x — > Ax. The superscript then stands for the power of 1/A). It is easy to 
show that k^\x) = p UJ (x), the classical momentum, solution of the Hamilton- Jacobi equation ([7]). 
The second equation is not less remarkable: it is a total derivative and it has the universal form 



:d x In 



uv(x) - v 2 p Ld + \d p F 2 (pu) 



where v gr = 1/d^p^ is the group velocity. Since k^ is purely imaginary, it governs the mode 



-d x In [F(p Ul) Vgr 



(A4) 



amplitude. Eq. (Al) constructed with k w = p, 

WKB 



few gives the generalized x-WKB expression 



(X) 



(A5) 



Its Fourier transform evaluated at the saddle point gives the p- WKB mode of Eq. ( 54 ) . 



Eq. ( A4 ) guarantees that the scalar product of Eq. (11) evaluated with Eq. ( A5 ) is exactly 



conserved. Eq. (A2) also guarantees that the development of k is alternated: even terms are real, 
while odd ones are imaginary. 



2. Leading deviations due to dispersion 



In this paper, we are interested in solving Eq. Q in the limit of weak dispersion, that is, 
for A large enough. The aim of this section is to precise the meaning of A large enough' by 



computing the scaling of the errors made when building the WKB basis of Sec II A 2 Two different 
approximations have been used. First we built approximate solutions of Eq. (8j) with the above 
WKB modes. Second, we solved Eq. ([7]) in the limit of large A in Sec II A 1 and we used the 
approximate roots to compute the WKB modes. 



To proceed, we estimate the next order term of Eq. (A3) in the limit of weak dispersion (WD), 



i.e. by dropping terms of order 2 in 1/A. In this regime, k^ is solution of 



2[F(p u )v gr {p u )] k$ D (x) 



i8 x + k^] ■ \k^d k (F Vgr ) 



(A6) 



Because we master the modes near the horizon in the p-representation (see Sec II B and III), we are 
only interested in the error accumulated from infinity, where the x-WKB approximation becomes 



exact, to ^pasting of Eq. (91), at the edge of the near horizon region. This error is estimated by 

(2) 

evaluating the integral of k\J D from x pas tmg till oo. Indeed the exact mode cj) u can be approximated 
by 



+ *(*)). 



(A7) 
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To evaluate <p^ KB we use Eq. <\A5h and the approximate roots of Sec 



II A 1 



errors governed by y of Eq. (38). Hence, near x pa sting the total error is 

/"^pasting , . 

er ^ / k\ v 1 D (x)dx +0(y(x pas t ing )). 



This introduces extra 



(A8) 



This error term behaves quite differently for in and out modes, hence we shall study it separately. 



For the in modes, solving the Eq. (A6), we get 

9v' 2 
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(A9) 



16A|l + f|2 4A|l + v|2 

Since this is not a total derivative, the integral depends on what happens all along the way 
from oo to £ pas tmg- However, when the profile v is smooth enough, the accumulated error is 
essentially 



A|l + u| 



+ 0(y p ) 



C pasting 
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For the out modes, the leading order correction arises from k^\ Indeed, in the limit A — > oo, 
the out modes are WKB exact because of conformal invariance and fei = 0. Therefore, 

A— >oo ' 

for finite A, k^ will be the dominant contribution to the error term of Eq. (A7). One finds 

6co 2 v' 
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This means that 
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Here the pasting happens on both sides, hence in the latter expression, one should understand 
Ap for and for y?° ut . As expected, the corrections to red-shifted out modes are 

subdominant with respect to those of in modes. 

It is interesting to notice that the validity of the x-WKB approximation in a given flow and 
for a given dispersion relation depends on the exact form of the mode equation associated with 
Eq. ([7]). Indeed, when the mode equation is not conformally invariant in the limit A — > oo there is 
an extra validity condition: 



< 1. 



(A13) 



It is due to the fact that for low frequencies, the left and right moving modes mix even in the 
absence of dispersion. This mixing was studied in [2S], and it was numerically shown that these 
effects stay (in general) subdominant. 
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Appendix B: Application: The black hole laser 
1. The setup 

When a flow contains both a black and a white horizon, and is subsonic on both asymptotic sides, 
HR is self- amplified in presence of superluminous dispersion [5HE2]) and this leads to a dynamical 
instability [53]. When the horizons are well separated, a powerful '5- matrix' description can be 



applied to describe the mode propagation. In this picture, as in Eq. (12), the matrix acts on the 
two right moving modes (tp u , (ip-. u )*). Because the horizons are well separated, it is legitimate to 
factorize S as S = SV^SbS^, where Sw and Sb are the 'on-shelF Bogoliubov matrices associated 
with the white and black hole respectively, and S2 and S4 contain phase shifts characterizing 
Hamilton-Jacobi propagations from the white hole horizon to the black one. These have the form 

*=' ^J' (B1) 

* (B2) 




where designates S u of Eq. (32) evaluated for the root k™, and similarly for the other two 



actions. The names of the roots are those associated with the black hole horizon, so that their 



expressions can be found in Sec II A 1 We notice that the precise values of the end points of 
integration in these three actions depend on the phase conventions that have been adopted to 
compute Sw and Sb, so that S is independent of these arbitrary choices. Notice also the minus 
sign before S4 because the corresponding propagation is backwards in x. For a more detailed 
analysis, we refer the reader to [53]. Here, we only recall the main points. 

For a real frequency uj > 0, the eigen-mode is asymptotically a plane wave of positive norm 
taken to be unity. To compute b w , the amplitude of the negative norm mode trapped between the 
two horizons, we require that the mode be single valued. This condition gives 

(B3) 

where e 1 ^ is the transmission coefficient, here a phase shift. In addition to the real spectrum, 
there exists a discrete set of pairs of complex frequency modes which encode the instability. To 
construct these modes, we look for frequencies A a = u a + iT a (where a labels a discrete set) such 
that the mode is asymptotically bounded (in fact square integrable). As demonstrated in [53j, the 
unstable modes (r > 0) are purely owi-going and thus satisfy 
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A a 




(B4) 
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This equation has a solution if and only if 

S , 22 (u; + «T) = 1. 

This is the algebraic equation that fixes the set of complex frequencies A a . 



(B5) 



Predictions 



To compute these complex frequencies, we perform an expansion in the (3 U coefficients in Sw 



and Sb- To zeroth-order in these coefficients, Eq. (B5) gives a Bohr-Sommerfeld condition that 
fixes ui a , the real part of the frequency: 



B 



{p°^(x) - p m Jx)) dx - arg(a^a^) = 2vrn, n 6 N*, (B6) 
where the end point x^ (respectively xf p ) refers to the turning point near the white hole (respec- 



tively black hole) horizon. The l.h.s. of Eq. (B6) displays both standard and unusual features. 



On the one hand, as expected, one finds the classical action evaluated along a closed loop from 
one turning point to the other one. One the other hand, one does not find the usual phase shift 
(= 7r) that accounts for the two reflections when dealing with Schrodinger type problems where 
the modes near the turning points can be well approximated by Airy functions. Indeed, using 



Eq. (76) and Eq. (71), one sees that — arg(a ( |j a B ) differs from tt. To characterize the difference, 



we introduce 

eH=arg(a^df)+vr, (B7) 

= arg(f (u/kw)) + arg(f (w/kb)). 

Using Eq. (73), one sees that it is only in the limit uj/k — > oo that one recovers the standard result, 
i.e. e = 0. For smaller values of lo/k, e accounts for the non-trivial phase shift due to the reflections 
on the two horizons. It arises from the fact that the reflected modes cannot be well approximated 
by Airy functions, something not discussed in [M] . 

In Fig. [7] we have compared the numerical results obtained using the code of [3S] with the 
theoretical predictions evaluated with and without e. The improvement of the agreement when 
including e is clear. 



To second order in Eq. (B5) fixes the growing rate T for a given value of uj: 
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cos(V>), 



(B8) 



where T u is the classical time for the trapped mode to make a round trip between the two horizons 

T u = d u (5° u j - Sil + n- e(w)) , (B9) 



The writing of Eq. ( B6 \ is obtained only if the Bogoliubov coefficients a™ and a" of Eq. ( 12 1 are computed with 



a phase 80 of Eq. (33 1 common to in and out modes. It is now clear why it is convenient to work with this choice. 
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Figure 6: Evolution of the real part (left plot) and the imaginary (right plot) of a complex frequency as a 
function of L, the distance between the 2 horizons. These curves have been obtained by making used of the 
numerical techniques described in |38j . The parameters used are K w /nb = 0.5, D = 0.5 and A/ Kb = 8, and 
we consider the n — 22 discrete mode. 



0.014 | 1 1 1 1 1 1 1 




Figure 7: Relative errors between the numerical results and our theoretical formulae of the real and imaginary 
parts of the same complex frequency as in the above Figure. The continuous lines take into account e of 



Eq. (B7|, while the dotted lines are based on the standard expression e = 0. The improvement of the 



estimation is clear, and therefore the necessity of computing the phases in Eq.(76l is established, 
and the phase ip reads 

^ = ^-^ + arg^-^J. (B10) 

This non-trivial phase arises from the interference between the two pair creation amplitudes, that 
occurring at the white hole and that at the black hole. A careful calculation gives a remarkably 
simple result 

( /-Pmax "I 

V> = Re / (X*(p) - X^{p)) dp\+7r. (Bll) 
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The constant phase shift follows from Eqs. (76 82) which give 



arg 



TT. 



(B12) 



The contribution from the classical actions is best expressed in p-space using the function X u {p) 



of Eq. (22). However, because the flow profile v is not monotonic in the black hole laser setup, 
X w is no longer unique. Assuming that v < has a single minimum at x = between the two 
horizons, we define X^f (p) (resp. X^{p)) as the solution of Eq. (22) of negative values describing 



the propagation towards the white hole (resp. positive values associated with the black hole). 
Both of these semi-classical trajectories run from a positive maximum value p max = (0) to a 
minimum negative value p m in = p™ w (0). In Eq. (Bll) we took the real part of the integral in 
order to remove the imaginary contributions (= ittu/kw and ittlo/kb) that arise when p flips sign, 



see the discussion after Eq. ( 33 ) . The contribution for p > accounts for the propagation of the 



positive norm mode, whereas that with p < for that of the trapped mode. What is remarkable 
is that when they are combined, the net result for 5™ — S^ u takes the form of the first term of 



Eq. (Bll). Notice again that this simple form is found only when using mode bases such that 



Eq. (B12) applies. 



Interestingly, Eq. (Bll) has the same structure as Bohr-Sommerfeld condition of Eq. (B6) with 



the role of x and p interchanged, i.e. Eq. (Bll) is a closed loop in p-space of some X UJ (p). What 
is unusual is that the action receives imaginary contributions when p changes sign. In fact when 
removing the restriction to the real part of the action, one gets 



3 B 
a w a B 



(B13) 



In other words the ratios /3 w /a w can be blamed on, and therefore computed from, the imaginary 



contributions of the action S u of Eq. (32) that arise when p flips sign. An early version of this 
relation was used for relativistic modes in [48J, and it is at the core of the so called 'Unruh' 
modes [39]. It was adapted to dispersive waves in [16] and implicitly used above when computing 



the connection coefficients in Sec. Ill It was also recently exploited in [50J in a similar context. 



To conclude, we underline that unlike the robustness of Hawking spectrum, Eq. (B6) and 



Eq. (Bll ) are not necessarily valid for very small frequencies. Indeed, we have made two additional 



approximations. First that the S- matrix can be factorized, i.e. that the propagation is correctly 
described by separating the scattering on each horizon and by WKB propagations between the 
horizons. Second, that the growing rate of the unstable modes are small (TT^ <C 1) when developing 
Eq. (B5). This condition is violated when uj — > because (3/ a becomes of order 1. However, when 



the density of unstable modes is not too low, the numerical analysis confirms that the above 
expressions correctly characterize the complex frequencies, as can be seen from Fig(6]and[7] 
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Appendix C: Massive case 



In this appendix we add a mass m on top of the ultraviolet dispersion of Eq. ^ : 



F 



m + ( p + 



2A 2 



(CI) 



and focus on the novel features brought in by the mass. We suppose that the two scales are well 
separated, i.e. m <C A. In this limit one has 



Fmijp) = V m2 +P 2 + 



P 



2A 2 y / m 2 + p 2 

Moreover, since the last term becomes non negligible only when p ~ A > m, we get 

^3 



(C2) 



F m (p) = Vm 2 + P 2 + ^. (C3) 
Therefore, to first order in m/A, F m is a sum of the relativistic massive dispersion plus a dispersive 



term. In the near horizon region, the mode in p-space is still given by Eq. (55), and, as in Sec. Ill 
the various modes in x-space are given by contour integrals 
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(C4) 



(C5) 



encodes the modification of the phase due to the mass. As before, the choice of the contour C 
dictates which mode one is considering. 



In the following, we construct the generalization of the decaying mode of Sec. IIIA1 because 
this is enough to extract the Bogoliubov coefficients. To this end, we must choose a branch cut to 



define both the ln(p) and \J m 2 + p 2 appearing in G m . These functions introduce three branching 
points, p = 0, and p = ±im. Here, we take the line — iR + extended until im to be the branch cut, 
as shown on Fig. |8j Hence in the limit m — > we recover what we did in the body of the paper, 
see Fig. [4| To compute Eq. (C4), we proceed as in Sec III A 1 



When x > 0, the introduction of a mass does not alter the discussion of Sec. IIIA1 since the 
saddle at p s = ih.y/2nx is well above the singularity at im for x sufficiently far away from the 
horizon. Hence 



JG m (p s ) 



x <PL x 



(C6) 



where ipZ) is the massless dispersive decaying mode of Eq. (46). Moreover, we have G(p s ) ~ 
because p s S> m in the region of interest. For x > 0, the mode is thus rapidly decaying, on a scale 

2 

kx ~ (re/A)s, and in the relativistic limit (i.e. A — > oo) it vanishes. Therefore, as in the massless 
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Figure 8: Representation of the contours in the p-plane to get (0°^)*, the out mode of negative norm. The 
hatched regions are the asymptotically forbidden ones. The wavy line is the branch cut of ln(p), and the 
bold line is what must be added because of the mass. The left panel is valid for x > 0, and the right one 
for x < 0. 



case, this mode is proportional to 



iout ■> 
— aj,m,/ 



*, the negative norm out going mode. Indeed, if it were 



containing a small amount of the positive norm out going mode 4 > 'umi ^ would oscillate on the 
right side of the horizon until kx ~ (co/m) 2 3> (k/A)5 which gives the location of the turning point 
where the mode is reflected due to its mass, or due to a perpendicular momentum [55| [56] . 

For x < 0, using the analytic properties of <fiuj(p), we deform the contour C into the union of 
C%, Ci and C3 shown in Fig [8j On C\ and C2, there are two saddle points at p = ±A^/ k\x\ that 



describe the high momentum incoming modes, as in Sec. IIIA1 Their contribution is 



c,uc 2 



urn • 3n , 

(e « e 1 4 



+ e L i x Lpl 



(C7) 



where are the massless dispersive in modes of Eq. ( 42 ) and Eq. ( 43 ) . Moreover, the saddle point 



approximation is controlled by the parameter A(x) of Eq. (34), irrespectively of the mass m<A. 
Hence, as expected, the high momentum contributions of the out mode are mass independent. 
Along C3 instead, we perform a strong limit A — > 00 as in Eq. (64), and we get 



(X) 



1 



v/4 



TTK JC 3 



P 



-l e i(xp+G m (p)) 



dp 



\/p 2 + m 2 



2tt* 



(C8) 



which is a massive relativistic mode of negative norm. Up to a complex amplitude A^, it gives 
(99°^; m )*, the low momentum out branch of the globally defined mode (</>°^ m )*. In brief, the 
mode obtained with the contour C is (c/>°^ m )*: For all w, it decays for x > 0, and on the left 
side, it contains three WKB branches 



A, 



^out \ 



6* 



(C9) 
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Since (c^™ w )* and cp^ are normalized and have opposite norms, their relative coefficient furnishes 
the ratio of the near horizon Bogoliubov coefficients \$u,m/okj,m\- 5 From Eq. (C7) and Eq. (C9), 



we obtain \Pu m /cx u m \ = e -7RJ//K . It is independent of m and has the standard relativistic value. 

We decided to study (4>°^ m )* m order to be able to discuss the asymptotic Bogoliubov coeffi- 
cients. So far indeed our calculation is restricted to the near horizon regions l.b of Fig. [2} Had we 
studied the positive norm mode ^ ut , we would have faced the propagation in the regions 2 and 3 
on the right of the horizon. Below the critical frequency discussed in footnote [5j the mode ^° ut 
is completely reflected. Above that threshold, (/>° ut is partially reflected in a non universal manner 
that depends on the actual profile v(x). However since this scattering is elastic, as it mixes left and 
right moving positive norm modes (these being the only ones present in that region), it does not 
affect the pair creation probabilities of asymptotic quanta, i.e. it fixes 'greybody factors'. In other 
words, the propagation in the right region has no influence on pair creation probabilities, both for 
to below and above to^. 

On the contrary the mode mixing on the left regions 2 and 3 will affect them because this mixing 
involves modes with norms of both signs. For instance, in Minkowski space-time no asymptotic 



particle is created even though Eq. ( C9 ) applies near the Rindler horizon described by Eq. (j8j) with 

v (x) = -(l-2/ex)3. (CIO) 



This implies that the far away scattering undoes the one occurring near the horizon, as can be 
verified by direct calculation, using Bessel functions. For the smooth and asymptotically constant 
profiles we consider in this paper, WKB modes nevertheless provide a controlled approximation on 
the entire left region. Indeed, the errors due to the high momentum modes are m independent and 



small whenever Ap" of Eq. (92) obeys Ap" 3> 1. The error due to the low momentum mode can be 



bounded using the treatment of AppjAj For the reflected modes, for ^ to < cj^, one gets 

p /oexact 



K 



0(-j. (Cil) 



/5cxact 

Because the error grows like 1/to for lo — > 0, there is no guarantee that the extra mode mixing will 



not significantly affect the values of the coefficients given in Eq. ( C9 ) . Hence further investigation 
is necessary to understand under which circumstances the standard divergence of |/3 W | 2 in n/to 
is recovered in the massive case. We are currently examining this interesting question. Let us 
conclude by noticing that since the asymptotic value of the wave vector k° ut is finite for to — > 0, 
this opens the possibility of having a phenomenon similar to what we found when studying white 



We have used the same conventions as in the massless case, see Eq. ( 12 1, because the relevant part of the scattering 
matrix is still 2x2, for all values of to. This is non trivial because the dimensionality of the asymptotic modes is 
different from that of the massless case. Indeed, when there is a mass, the left going mode 0jj ft only exists above a 
critical frequency oj^ = mil — {v^) 2 ) 1 ^ 2 where v^ s is the asymptotic velocity on the right side. However, as later 
discussed in this Appendix, the mode mixing with </>Jj ft does not affect the asymptotic Bogoliubov coefficients. 
Hence <^J ft remains a spectator mode. 
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